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Vacuum polarization and photon mass in inflation 
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We give a pedagogical review of a mechanism through which long wave length photons can become 
massive during inflation. Our account begins with a discussion of the period of exponentially rapid 
expansion known as inflation. We next describe how, when the universe is not expanding, quantum 
fluctuations in charged particle fields cause even empty space to behave as a polarizable medium. 
This is the routinely observed phenomenon of vacuum polarization. We show that the quantum 
fluctuations of low mass, scalar fields are enormously amplified during inflation. If one of these 
fields is charged, the vacuum polarization effect of flat space is strengthened to the point that 
long wave length photons acquire mass. Our result for this mass is shown to agree with a simple 
model in which the massive photon electrodynamics of Proca emerges from applying the Hartree 
approximation to scalar quantum electrodynamics during inflation. A huge photon mass is not 
measured today because the original phase of inflation ended when the universe was only a tiny 
fraction of a second old. However, the zero-point energy left over from the epoch of large photon 
mass may have persisted during the post-inflationary universe as very weak, but cosmological-scale, 
magnetic fields. It has been suggested that these small seed fields were amplified by a dynamo 
mechanism to produce the micro-Gauss magnetic fields observed in galaxies and galactic clusters. 

PACS numbers: 98.80.Cq, 98.80.Hw, 04.62. +v 



I. EXPANDING UNIVERSE AND INFLATION 

The universe is expanding, but with a rate so tiny that 
it can only be seen by spectroscopic analysis of stars in 
distant galaxies. Suppose the light from such a star con- 
tains a distinctive absorption line measured at the wave 
length A. If the same line occurs at wave length As on 
Earth, we say the star's redshift is, 



10 -3 the speed of light c. This motion gives rise to 
a special relativistic Doppler shift of Az ~ ±10~ 3 . If 
spacetime was not expanding, this shift would be the only 
source of nonzero z, and averaging over many stars at the 
same luminosity distance would give zero redshift. That 
is just what happens for stars within our galaxy. How- 
ever, the luminosity distances of stars in distant galaxies 
are observed to grow approximately linearly with their 
redshifts, 



A 



- 1 . 



(1) 



qo)z 



0(z 3 



(3) 



One can also measure the star's flux T . If we understand 
the star enough to know it should emit radiation at lumi- 
nosity C, we can infer its luminosity distance dr,, which 
is the star's distance in Euclidean geometry, 



T = 



1-Kd\ 



dr = 



c 



(2) 



Astronomers measure the expansion of the universe by 
plotting z versus dj, for many stars. 

Stars throughout the universe move with respect to 
their local environments at typical velocities of about 



'Electronic address: T.Prokopec@thphys.uni-heidelberg.de, Tomislav 
T Electronic address: woodard@phys.ufl.edu 



Equation J3J means that more and more distant objects 
seem to recede from us with greater and greater speed. 
A common analogy is to the way fixed spots move apart 
on the surface of a balloon that is being blown up. 

The constant, H ~ 2.3 x 1CT 18 Hz, is called the 
Hubble parameter. (When expressed in units used by 
astronomers, Hq — 71 km/s/Mpc, where 1 Mpc corre- 
sponds to the distance light traverses in about 3.26 mil- 
lion years.) Its name honors Edwin Hubble, who estab- 
lished the (nearly) linear relation Q] in 1929 based on his 
observations, and on earlier work of Slipher and Wirtz 0] . 
The other constant in Eq. © is known as the decelera- 
tion parameter, qq. Observations of Type la supernovae 
(whose luminosities can be precisely inferred) up to the 
enormous redshift of 1.7 indicate qo ~ —0.6 
Prokip^ciie'eiftieVT °f spacetime is described by a symmetric 
tensor held g^ v {x) known as the metric. It is used to 
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translate the coordinate labels of points = (ct, x) into 
physical distances and angles. For example, the square 
of the distance between x M and an infinitesimally close 
point x^ + dx^ is known as the invariant interval, 



ds 2 



g fjLV (x)dx 11 dx v 



(4) 



Note that we employ the Einstein summation convention 
in which repeated indices are summed over 0, 1, 2, 3. 

The transition from nearby stars, whose redshifts are 
dominated by local motions, to more distant stars which 
obey Eq. © , is known as entering the Hubble flow. It is 
typical in cosmology to ignore local features and model a 
simplified universe that has only the overall expansion ef- 
fect. Such a universe does not change in moving between 
spatial points at the same time, nor are there any special 
directions. The first property is known as homogeneity; 
the second is isotropy. 

With a simplifying assumption — about which more 
later — the invariant interval of a homogeneous and 
isotropic universe can be written as 



c 2 dt 2 



a 2 (t)dx ■ dx . 



(5) 



From this relation we see that t measures physical time 
the same way as in the Minkowski geometry. However, 
the spatial 3-vector x must be multiplied by a(t) to give 
physical distances. For this reason a(t) is known as the 
scale factor. Its time variation gives the instantaneous 
values of the Hubble and deceleration parameters, 

. i Oj , . ad H 

H(t) = ~, q(t ) = - — = -!- — , 6 

The subscripts on Hq and qo in Eq. and the sub- 
sequent discussion indicate the current values of these 
parameters. 

Homogeneity and isotropy restrict the stress-energy 
tensor to only an energy density p(t) and a pressure p(t) , 



T, 



(K) 



-p(t)goo, T(h = 0, Tij = p(t)g ij , (7) 



where i and j are spatial indices. In this geometry Ein- 
stein's equations take the form, 



3H 2 = 8irGc~ 2 p, 
-2H-3H 2 = 8ttGc~ 2 p, 



(8) 
(9) 



where G is Newton's constant. The current energy den- 
sity is, 



p a = ^ ~ 8.5 x 10- 10 J/m 3 



(10) 



equivalent to about 5.7 Hydrogen atoms per cubic meter. 
If we solve for the instantaneous deceleration parameter, 



1 3p 
9= 2 + V 



(11) 



we find that po — — 0.7/?o 



By differentiating Eq. JSJ) and then adding 3H times 
Eq. (JSJ) plus Eq. we derive a relation between the 
energy density and pressure known as stress-energy con- 
servation, 



p= -3H(p+p). 



(12) 



If we also assume a constant equation of state, w = 
p(t)/p(t), Eq. itT^I) can be used to express the energy 
density in terms of the scale factor, 



P(*)=px(^) 



a(t)\-3(i+™) 



(13) 



The substitution of Eq. O in Eq. $ gives an equation 
whose solution is, 



a(t) = a\ 



l + -(l + u;)JTi(t-ii) 



(14) 



The cases of w = +-|, 0,— and —1 correspond to radi- 
ation, non-relativistic matter, spatial curvature, and vac- 
uum energy, respectively. (We do not discuss here scalar 
field matter 0, usually referred to as quintessence 0, 
with nonconstant w.) The cosmology for each pure type 
of stress-energy can be read off from Eqs. <|13[) and 114|l . 



Radiation 
Non-Relativistic Matter 
Curvature 
Vacuum Energy 



p oc a 
p cx a~ 
p oc a~ 
p cx 1, 



a(t) cx (Hit)*, 
(15) 

a(t) cx (Hit)i, 
(16) 

a(i) cx Hit, 

(17) 

a(t) cx e Hlt . 

(18) 



The actual universe seems to be composed of at least 
three of the pure types, so the scale factor does not have 
a simple time dependence. However, as long as each type 
is separately conserved, we can use Eq. (|13H to conclude 
that 



P(t) 



Prad Pmat 



a 4 (<) a 3 (i) a 2 (t) 



Pv 



(19) 



As the universe expands, the relative importance of the 
four types changes. Whenever a single type predomi- 
nates, we can infer a(t) from Eq. 1)14(1. This different de- 
pendence is one reason it makes sense to think of an early 
universe dominated by radiation, Eq. I|15l) . evolving to a 
universe dominated by non-relativistic matter, Eq. 116( 1 . 
It is also how we can understand that the current uni- 
verse seems to be making the transition to domination 
by vacuum energy, Eq. 1(18(1. 

Under certain conditions there can be significant en- 
ergy flows between three of the pure types of stress- 
energy. For example, as the early universe cooled, mas- 
sive particles changed from behaving like radiation to 
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behaving like non-relativistic matter. This change would 
increase p mat and decrease p ra( j in Eq. (|19|) . The param- 
eter that cannot change is that of the spatial curvature, 
p cur . Strictly speaking, we should not regard spatial cur- 
vature as a type of stress-energy, but rather as an addi- 
tional parameter in the homogeneous and isotropic met- 
ric Eq. JSJ. We have avoided this complication because 
the extra terms it gives in the Einstein equations ©-© 
can be subsumed into the energy density and pressure, as 
we have done, and because the measured value of p C ur/ a o 
is consistent with zero 0, @ . 

The cosmology in which a radiation dominated uni- 
verse evolves to matter domination is a feature of what 
is known as the Big Bang scenario. Although strongly 
supported by observation, the composition of p at the 
start of radiation domination (t — t r and a = a r ) does 
not seem natural, 

Prad a r Pvac ^ Pcur 

a r . (20) 

It might be expected instead that each of the three terms 
was comparable, in which case the universe would quickly 
become dominated by vacuum energy. There is no ac- 
cepted explanation for the first inequality in Eq. (|20l) or 
for the seeming coincidence that p ma taQ 3 ~ p vac . How- 
ever, the second inequality in Eq. (|2(Jfl finds a natural 
explanation in the context of inflation. 

In 1980 Alan Guth @ suggested that the Big Bang sce- 
nario was preceded by a period of vacuum energy dom- 
ination, or inflation, following which the vacuum energy 
changed almost completely into radiation. (Cosmologies 
that include a period of vacuum energy domination were 
independently considered by Starobinsky 01, Sato [TT| . 
and by Kazanas 0.) Suppose that all types of stress- 
energy are equally represented at some very early time. 
We see from Eq. (I19fl that the total energy density rapidly 
becomes dominated by vacuum energy, following which 
the scale factor grows exponentially with a constant Hub- 
ble parameter, Hj. 

The duration of inflation in units of 1/ffj is known as 
the number of inflationary e-foldings Nj. Viable models 
must have Nj > 50, and much larger values are common. 
If Pcur/ct/ ~ Pvac at the start of inflation, Eq. l|19fl shows 
that the curvature is negligible at the end, 

P^IA^ C^ 2 =e -2JVz £ 10 -44 (21) 
Pvac ^ a r ' 

Inflation makes the other types of stress-energy even 
smaller, but there are mechanisms through which vac- 
uum energy can be converted into radiation. This pro- 
cess, which we will not discuss, is known as reheating. 

Inflation also explains how the large scale universe be- 
came so nearly homogeneous and isotropic. This expla- 
nation is crucial because gravity makes even tiny inho- 
mogeneities grow, and the process has had 13.7 billion 
years to operate. It is believed that the galaxies of to- 
day's universe had their origins in quantum fluctuations 
of magnitude Ap/p ~ 10 -5 , which occurred during the 



last 60 e-foldings of inflation. The imprint of these fluctu- 
ations in the cosmic microwave background has recently 
been imaged with unprecedented accuracy by the WMAP 
satellite H®. 

The fact that WMAP did not see the imprint of quan- 
tum fluctuations of the metric field sets an upper limit 
of Hi < 3.4 x 10 38 Hz. No one knows what caused in- 
flation, but a common assumption is that it occurred at 
the grand unified energy scale Sgut — 10 16 GeV ~ 10 6 J 
at which electromagnetic, weak, and strong interactions 
attain equal strength. From Eq. JSJ) this implies, 

where we used pgut = -Egut/Vgut, Vgut = 
(Tlc/Egjjt) 3 defines the GUT energy scale volume, and 
h ~ 1.05 x 10~ 34 Js is the reduced Planck constant. 



II. VACUUM POLARIZATION IN FLAT SPACE 

Flat space corresponds to a(t) = 1 in Eq. 

dsflat = rj^dx^dx 1 ' = —c 2 dt 2 + dx ■ dx . (23) 

Note that the zero component of a spacetime point x^ = 
(ct,x) is x° = ct, so all components of = d/dx 11 have 
the dimension of inverse length. Repeated Greek indices 
are summed over 0,1,2,3 — for example, d 2 = d^d^ - 
whereas repeated Latin indices are summed over 1,2,3 
- for example, V 2 = did L . A dot denotes contraction 
over the appropriate index set, for example, k ■ x = fc M cc M 
and k ■ x ^= k^x^. 

To make Lorentz invariance manifest, we employ the 
Heaviside-Lorentz system of units in which Maxwell's 
equations take the form, 

V-E = p, V x B- d E = c' 1 !, (24) 
V-B = 0, VxE + d o B = 0. (25) 

Here E(t, x) and B(t, x) denote the electric and magnetic 
fields, and the charge and current densities are p(t,x) 
(for this section only) and J(t,x). The more familiar 
MKSA formulation of electrodynamics follows from the 
substitutions, 

E — > eo-Emksa, B — > » /— B UK3A , and c — > , 
V Po v e °^° 

(26) 

where eo and po are, respectively, the electric permittiv- 
ity and the magnetic permeability of free space. 

It is well known that Eq. (|25|l can be enforced by repre- 
senting the fields using a vector potential = (A , A.A, 

E l = d A, - diA , B i = -e^ k djA k . (27) 

Equations l|24|) combine to the relativistic form, 

d u F u i* = c^J" , (28) 
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using the field strength tensor = d^A u — d„A^ and 
the current 4- vector = (cp, J). 

Material media such as air and glass consist of an 
enormous number of atoms with negatively charged elec- 
trons bound to positively charged nuclei. On macroscopic 
scales such a medium appears neutral and free of cur- 
rents, but the application of external fields can distort 
the bound charges to induce a density of atomic electric 
dipole moments known as the polarization P(t, x), which 
we illustrate on a gas of polarized atoms in figure ^ Av- 
eraging the actual charge density to remove its violent 
fluctuations on the atomic scale leaves whatever charges 
are free, minus the gradient of P, 



(p) = Pit, 



V-P. 



(29) 



The medium's density of atomic magnetic dipole mo- 
ments is known as its magnetization M(t,x). A similar 
averaging of the current density gives, 



(J) = J frcc + cd P + cV x M . 



(30) 



Moving the polarization and magnetization terms to the 
left-side of Eq. l|2"4"|l leads to the macroscopic Maxwell 
equations, 



V • D = jOfrec 
V X H~d D = C^Jfre 



(31) 



where D = E + P and H = B - M. 

Linear, isotropic media with no frequency or wave 
number dependence are characterized by, 



P = XeE, 



M = 



Xr, 



1 + Xr, 



-B. 



(32) 



The dimensionless parameters Xe and Xm are known 
as the electric and magnetic susceptibilities. We would 
like to express Eq. 1)31(1 as a single tensor equation like 
Eq. H28|) . For the case of constant susceptibilities the 
result is simple: 

d v F v » + W u 'A v =c- 1 Jg„. (33) 

where IF"" is the following tensor differential operator, 

W v = Xe(v^d 2 - d^d") 

- (Xe + T §z r )rrv vj (%V 2 - didj) . (34) 



The polarization tensor encapsulates the medium's ef- 
fect on electromagnetic forces and on propagating elec- 
tromagnetic fields. It is useful to recall the familiar for- 
mulae for the relative permittivity and permeability and 
for the index of refraction, 



e = 1 



Xe 



P = l+ Xr. 



n = (efj,)- 



(35) 



The susceptibility Xe in Eq. I|34|) gives the medium's cor- 
rections to the electric response to a static distribution of 



charge. Positive Xe means that the medium's dipoles line 
up to weaken an applied electric field by a factor of 1/e. 
This effect is known as charge screening. The other term 
in Eq. I|34|l can be understood by recasting its integrand, 



Xe 



1 + Xr, 



1 

e 



l)-=Sn 2 
P 



(36) 



Together with e, Sn 2 gives the medium's corrections to 
the magnetic response to currents. It also governs the 
speed c/n at which electromagnetic waves propagate, 
such that, whenever Sn 2 ^ (n ^ 1), the speed of light 
differs from that in the (Minkowski) vacuum. 



i 



X 



FIG. 1: A gas of polarized atoms. In the absence of an exter- 
nal electric field the dipoles orient randomly (FIG.^i.). When 
an external field E is applied, the dipoles tend to line up with 
it (FIG. ^i). This alignment produces a net polarization, 
P = XeE> which weakens the electric force by 1/(1 + Xe). 
There is a similar effect in vacuum due to oppositely charged 
pairs of evanescent, virtual particles. 

The susceptibilities of real media typically vary accord- 
ing to the frequency and sometimes even the wave num- 
ber of the external field ^| . One reason for this variation 
is that the bound charges in a medium cannot respond in- 
finitely quickly to changes in the applied electromagnetic 
fields, so the polarization must be attenuated at high 
frequencies. There also can be resonant effects when the 
applied electromagnetic fields perturb systems of bound 
charges near their natural vibrational frequencies. 

Whatever the reason for it, we should understand that 
frequency and wave number dependence invalidates the 
local constitutive relations l|32|l and hence, the tensor 
equation l|33l) we have derived from them. Maxwell's 
equations in media (|31|l are still correct, but the rela- 
tions between the polarization and magnetization and 
the applied fields are only multiplicative in frequency- 
wave number space. To calculate P(t,x) and M(t,x), 
we first decompose the electric and magnetic fields into 
their amplitudes for each wave 4- vector = (w/c, k). 
This decomposition is equivalent to taking the Fourier 
transform, and we denote the result with a tilde as fol- 
lows, 



E\uj,k) 



cfixe™ 1 - 



//»••;•' 



E\t,x). 



(37) 



To get the polarization, we multiply by the k^ dependent 



5 



susceptibility and use the Fourier inversion theorem, 
d 4 k 



P l (t, 0)=J „e*»xe(*)^(w, k) . (38) 

The analogous procedure gives the magnetization. 

We have just seen that properly accounting for the 
dependence of real media requires us to perform two 4- 
dimensional integrations: first over spacetime to Fourier 
transform the fields, and then over This is tedious 
and unnecessary. By calculating the Fourier inverse of 
the susceptibility once and for all, 



Xe {x, X ) — 



d 4 k 



Xe(k)e 



ik-(x-x') 



(39) 



we can reduce the process of dealing with different ap- 
plied electric fields to a single spacetime integral, 



P l (t,x) = J d 4 x' X e(x,x')E l (t',x') . (40) 

An analogous simplification can be made for the magne- 
tization. 

We are now ready to give the generalization of Eq. ilMMI) 
which applies for the case of fc M dependent media, 



r frfd p F m {x)+ J d 4 x'[»Wlx, x')A u {x')= c~ l J»Jx). 

( 41 ) 

The polarization bi-tensor has the general form, 

+ rf i rf*[8 ii ff k d h -& i d i ]&n i {x, a !), (42) 
where <9 M — d/dx* 1 , d' = d/dx ,fl and, 



8n 2 (x, x') 



d 4 k ( Xm(fc) , 



ik-(x — x ) 



(43) 



We have actually done a little better than intended. For 
although Eqs. (|39|l and (|43(1 apply only to translation 
invariant and infinite media, Eqs. (|41|l (|42|l give the re- 
sponse from any linear, isotropic medium. For exam- 
ple, we could use this formalism to describe a system in 
which the local index of refraction varies in space or even 
in time. Although we shall not consider spatial depen- 
dence, this natural way of incorporating time dependence 
will be quite useful when we consider electrodynamics in 
an expanding universe. 

The designation of [TP] (a;, x') as a bi-tensor derives 
from general relativity in which the index /i transforms 
according to the vector space at x^ and the index v ac- 
cording to the vector space at x' 1 *. Note further that 
[TP] is transverse on both indices, 



d^U»}{x,x') = = d' v [^]{x 1 x') 



(44) 



Let us explore the dynamical consequences of an infi- 
nite, translational invariant medium for which Eqs. I|39|) 



and (|43() pertain. In this case we may as well Fourier 
transform (|41|l in the Coulomb gauge, kiAi(k) = 0. The 
fx = component, 



e(k)\\k\\ 2 A°(k) =c- 1 J°(k), 



(45) 



determines the scalar potential from the charge density. 
As claimed, the medium screens the electric forces by a 
factor of \/e{k). The fi — i equations are more interest- 
ing: 



e(k) 



2 (k) 



A\k) 



J j (k). (46) 



In addition to the response to a current, the 3- vector 
potential can also support plane waves that obey the fol- 
lowing dispersion relation, 



e(k) n- 2 (fc)||fc|| 



-2 2 
C UJ 



0. 



(47) 



Einstein's great contribution to quantum theory was the 
inference (from the photoelectric effect) that light is 
quantized in discrete photons of energy E = huj and 3- 
momentum p = Hk. 

When e(k) is nonsingular, Eq. (|4*7|) implies the usual 
relation, w — (c/n)\\k\\. For this case the energy vanishes 
as the wave length becomes infinite. However, suppose 
the medium obeys, 



Xe(fc) 



ra^c 2 
% 2 k-k' 



n{k) 



1 



(48) 



where, as we will see in the following, m 7 denotes the 
photon mass. Although the full significance of the sin- 
gular behavior l|48|l will become clear later, here we note 
that it may arise, for example, in the presence of charged 
particles that are constantly created, and which propa- 
gate with the speed of light. From the point of view of 
the photons, a large number of these particles may lie on 
its past light cone, at which k ■ k = 0, and may thus in- 
duce a large (singular) contribution to the susceptibility. 
The substitution of Eq. gHJ in Eq. g7J| gives, 



e(fc) n- 2 (fc)||fc|| 



ur 



= k-k 



2 2 



= 0. 



(49) 



Such a photon's energy is that of a massive particle, 



E = hu = 



IpII 



m 2 c 4 . 



(50) 



We have so far discussed classical media. Quantum 
field theory predicts that particlc-antiparticlc pairs are 
continually being created. They live for a brief period 
of time and then annihilate one another. The lifetime 
of such a virtual particle pair is governed by its energy 
through the energy-time uncertainty principle, 



At AE > h . 



(51) 
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The meaning of Eq. (|51|l is that a minimum time At is 
needed to resolve the energy with accuracy AE. Sup- 
pose each partner of a virtual particle pair has energy E. 
Before they emerged from the vacuum, the energy was 
zero, whereas it is IE afterward. This is an example of 
nonconservation of energy! However, Eq. 1)51(1 says that 
the violation is not detectable in a period shorter than 
H/2E, so virtual particles can survive roughly that long. 

All types of particles experience virtual particle cre- 
ation with all possible energies and 3-momenta. One 
way of understanding the electrostatic force is through 
the exchange of virtual photons. Because normal pho- 
tons are massless, they can have arbitrarily small ener- 
gies and can therefore survive long enough to mediate 
the force between distant charges. However, the life- 
times of massive particles are extremely short. For ex- 
ample, the minimum energy an electron can have is that 
of its rest mass, m e c 2 ~ 8.2 x 10~ 14 J. From Eq. ifST]) 
we see that an electron-positron pair can only live about 
h/2m e c 2 ~ 6.4 x 10 -22 s. Even moving at nearly the 
speed of light (which implies higher energy, and hence 
shorter lifetime), a virtual electron would only cover 
about 10 _13 m before annihilating, which is a thousand 
times smaller than the scale upon which the discrete 
electrons and nuclei are separated in atoms. This ex- 
plains why macroscopic experiments do not detect virtual 
electron-positron pairs. 

Charged virtual particles behave much like the bound 
charges of atoms in a polarizable medium. When no ex- 
ternal field is present, the positive partner of a virtual 
pair emerges as often in one direction as any other. How- 
ever, the application of an electric field makes it prefer- 
able for the positive partner to emerge in the direction of 
the field, while the negative partner emerges in the oppo- 
site direction. In this way even empty space can acquire 
a polarization. The effect is known as vacuum polariza- 
tion, and it is described with the same Eqs. (|41|l H42|) 
which we introduced to quantify the polarization of ma- 
terial media. 

Although all charged virtual particles contribute to 
vacuum polarization, the largest effect comes from the 
lightest particles because they live the longest. Elec- 
trons and positrons are about 200 times lighter than 
the next lightest charged particles, muons, so almost all 
vacuum polarization comes from them. By making use 
of rather sophisticated techniques of quantum electrody- 
namics, one can show that, at lowest order in a, they 
induce the following electric susceptibility [l4| . 



8wa I dxx(l 
'o 



d 3 p 

(27T) 5 



(52) 



where a = e 2 /Anhc ~ 1/137 is the fine structure con- 
stant. The integral over p represents the contribution 
from an electromagnetic field of wave vector k exciting 
an electron with wave vector p + xk and a positron with 



wave vector — p + (1 — x)k. 

The process described by Eq. 1(52(1 conserves 3- 
momentum (=fix wave vector), but it need not conserve 
energy, 



hck° — ► y m 2 c 4 + H 2 c 2 \\p + xk\\ 2 

+ \fm 2 e c 4 + h 2 c 2 \\p- (1 - x)k\\ 2 . (53) 

The energy of the electron-positron pair can be very 
much larger than the energy of the initial and final pho- 
tons, which means the pair can exist only a short time. 
As one might expect, pairs with very high \\p\\ do not con- 
tribute much susceptibility because they exist too briefly 
to polarize much. However, there are so many states with 
large \\p\\ that their net contribution diverges. That is 
why we have cut the integral off in Eq. I(52|) for \\p\\ < A. 
If we integrate over the momenta, and expand in powers 
of 1/A, we obtain, 



Xe(fc,A) 



4a 
n 
1 



dxx(l 



, /2hA 
In I 



) 



In 1 



( 1 + -( 1 -^S^)+°( 1 /A)]. (54) 



m e c/ 
k- 



We see that the susceptibility diverges logarithmically in 
the limit that A — > oo. This is a classic example of an 
ultraviolet divergence in quantum field theory. (The ad- 
jective ultraviolet is used because the problem originates 
from high energy — or ultraviolet — pairs.) 

The infinite susceptibility is not directly observable be- 
cause it is a constant, independent of the wave 4-vector 
fc M of the applied electrodynamic field. All observable 
quantities depend as well upon the equally constant, 
bare charges of particles whose motions reveal the field. 
For example, consider the force between two particles of 
charge ebarc that are held a fixed distance r from one an- 
other. Because the system is constant in time, we have 
k° = 0. Because the physical dimension of the system 
is r, the maximum response is for wave vectors of about 
||fc|| ^2ir/r. As the separation becomes larger, the signif- 
icant wave vectors become ever closer to zero. For very 
large separations the force is therefore 1/r 2 times the 
quantity, e^ are /47r[l + Xe(0, A)]. The case of large sepa- 
rations (on the scale of 10~ 13 m) is one that we can access 
quite easily, and most everyone who takes introductory 
physics performs such a measurement. Because the re- 
sult is finite, it follows that the ratio e 2 aro /[l + Xe(0, A)] 
must be a finite number that we call e 2 , the square of 
the measured charge. In other words, the divergence in 
the unobservable quantity e 2 aro cancels the divergence in 
the equally unobservable quantity Xe(0, A) so that their 
ratio agrees with what we measure. To leading order in 
a, the finite fc M dependent susceptibility, which remains 
when measured charges are used, is, 



7 



Xe (fc) = lim \ Xe (k,A)-Xe(0,A)} 



2a 



(55) 



— / dx x(l — x) In 1 + x(l — x) 



Ti 2 k-k 



(56) 



This discussion illustrates how the process of renormal- 
ization works in quantum field theory. We mention it 
only to explain why the finite remainder Eq. (|56J) can 
make electromagnetic forces stronger at short distances. 
Recall that the least energetic electron-positron pairs 
can only survive long enough to travel about 10~ 13 m. 
More energetic virtual pairs are limited to even shorter 
distances. This means that charged particles separated 
by more than about 10 -13 m feel the polarizations con- 
tributed by virtual pairs of all 3-momenta. But at separa- 
tions of less than 10 -13 m, the lower energy virtual pairs 
leave the electric field between the two charges before 
fully polarizing. The net effect is less charge screening 
than at large distances, and hence a relative enhance- 
ment of the electromagnetic force at short distances. 

This effect is known as running of the force law, and 
it is seen routinely in precision measurements. In high 
energy accelerators such as LEP at CERN and SLC at 
Stanford University, charged particles have been brought 
as close as ~ 10~ 18 m. The substitution of k ■ k = 
(27r/lCT 18 m) 2 in Eq. J2J gives X e - -0.023, or a 2.3% 
enhancement of the electromagnetic force. 

If the electron mass had been zero, we would see the 
electromagnetic force law run even at macroscopic dis- 
tances. In that case the renormalized e 2 would be 4ttR 2 
times the measured force at some R, and this length 
would enter the formula for Xe(k), 



Xe(k) 



m e =0 



3tt 



ln(/i 2 £;•£;), where fi 



R 
2tt' 



(57) 



For r > R we would measure the force to be smaller than 
e 2 /47rr 2 , whereas it would be greater for r < R. The 
experiment could be done using the apparatus depicted 
in Fig. |2 [H G]| . 

We have so far discussed only the term Xe(x,x'). It 
turns out that Sn 2 (x,x') is zero for all relativistic quan- 
tum field theories in flat space. This must be so because 
the tensor coefficient of Sn 2 (x, x') in Eq. (|42|l breaks the 
Lorentz symmetry between space and time. One conse- 
quence is that the index of refraction is one, so vacuum 
polarization does not modify the speed of light. Because 
the invariant element of an expanding universe Eq. (J5J 
also distinguishes between space and time, we might ex- 
pect that Sn 2 (x,x') ^ when the Hubble parameter is 
nonzero, and we will sec that this is the case. 

Because Eq. l|56*)) has no pole at k-k = 0, the vacuum 
polarization from quantum electrodynamics preserves the 
photon's zero mass. This turns out to be a slightly 
mass and dimension-dependent statement. In 1962 Ju- 
lian Schwinger showed that zero mass electrons in one 




FIG. 2: A sketch of the Maxwell-Cavendish experiment. Two 
metal concentric spheres of radii a and b are first grounded, 
then the outer sphere is raised to a high potential. If Cou- 
lomb's law were violated, the voltmeter V would show a non- 
zero voltage. 



spatial dimension would induce the following electric sus- 
ceptibility 



Xe(k) 



4a 
~k~k 



(58) 



(In two spacetime dimensions e 2 has the dimension of 
energy /length, which means that a has the dimension of 
length -2 .) A comparison with Eq. (|48|l implies a photon 
mass of m 1 = 2-Jah/c. We will see later that the ex- 
pansion of spacetime can also induce a nonzero photon 
mass. 

Electrons and positrons, both of which are spin 1/2 
particles, are not the only kinds of charged particles. To 
obtain the susceptibility contributed by other kinds of 
spin i particles, we simply replace m e in Eq. (|56*|l by 
the appropriate mass. Charged particles with spin zero 
- which are known as scalars — also entail replacing 
the factor of x(l — x) which multiplies the logarithm in 
Eq. (|56|l by (1 — 2x) 2 /8. The susceptibility of a zero mass 
scalar would be j times that of Eq. (|57|1 . It is actually 
simpler to express the polarization of massless particles 
in position space by performing the Fourier transform of 
the electric susceptibility, as indicated in Eq. (|39J) . The 
result for a massless, charged scalar is 20, 21], 



a 



967T 



;d 4 {9(At)6(AT 2 )[l - ln(^ 2 Ar 2 )]}, 



Xe(x,x') 

'■IU7T- i. 

(59) 

where At = t - t', At 2 = (t - t') 2 - c~ 2 ||£ - x'\\ 2 and 
v = c/R is the frequency scale of renormalization. Note 
that Eq. I|59l) is zero whenever the point x'^ lies outside 
the past light-cone of x^. This feature, which is a fun- 
damental requirement on any % e (a:,a;') and 8n 2 (x, x'), is 
known as causality. 
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III. VIRTUAL PARTICLES WITH EXPANSION 

Leonard Parker was the first to give a quantitative as- 
sessment of how the universe's expansion can affect vir- 
tual particles (lsf. The mechanism is that the partners 
of a virtual pair must cover more distance getting back 
together than they did moving apart. This causes them 
to stay apart longer. Under certain conditions they can 
become trapped in the Hubble flow and become pulled 
apart, leading to physical particle creation. The purpose 
of this section is to explain why the effect is strongest dur- 
ing inflation and for massless scalars which possess the 
special property of minimal coupling to gravity, about 
which more later. 

First consider how the energy-time uncertainty princi- 
ple Eq. H51|) generalizes to the homogeneous and isotropic 
geometry in Eq. JSJ. Just like photons, a general quan- 
tum mechanical particle is characterized by its wave vec- 
tor k, which points in the particle's direction of propa- 
gation and has magnitude 2n divided by the particle's 
wave length. Now recall from Eq. © that the physical 
length between two fixed spatial points is a(t) times their 
coordinate separation, ft follows that the physical wave 
vector is k p ^ = k/a(t). The 3-momentum of a quantum 
mechanical particle is Ti times its physical wave vector. 
Hence the energy of a particle with mass m and coordi- 
nate wave vector k is, 



E{t,k) 



h 2 c 2 \ 



k\\ 2 /a 2 {t). 



(60) 



This changes with t so the energy-time uncertainty prin- 
ciple says we cannot detect a violation of energy conser- 
vation at time t + At from a pair of such particles created 
at t, provided that 



t+At 



dt' 2E(t', k)<h. 



(61) 



We see from Eq. H60JI that the growth of a(t) always 
reduces the energy relative to the constant scale factor. 
From Eq. 1|6I[) we see that the growth of a(t) always in- 
creases the time a virtual pair can survive. For a given k 
and time dependence a(t), the rate at which E(t, k) falls 
increases as the mass decreases. Hence, massless virtual 
particles experience the largest increase in their lifetimes. 

To understand why inflation maximizes the effect, con- 
sider the form of Eq. (|6f|l for a massless particle, 



2frc||fc|| 



t+At 



dt'- 



1 



o(f) 



<h. 



(62) 



For the radiation dominated scale factor Eq. I|f 5fl the in- 
tegral grows like (At) 2 ; for matter domination, Eq. (|f 6fl , 
its growth is like (At) 3; and the growth is logarithmic 
for curvature domination Eq. (| 1 7|) . In each of these cases 
the inequality is eventually violated as At grows. How- 
ever, for the inflationary scale factor in Eq. I|18|) . the inte- 
gral approaches a constant as At becomes infinite. This 



means that a long enough wave length pair need never 
recombine. Models of inflation are typically well approxi- 
mated by a(t) oc e Hlt , for which the bound Eq. Ij62|) takes 
the form, 



I All 



a(t) 



e 



— Hi At 



(63) 



Therefore massless particles of coordinate wave vector 
k are created during inflation whenever a virtual pair 
emerges with c||fc p h| ~ Hi. This condition on the wave 
number is known as first horizon crossing. 



inflation 


radiation era 


matter era / 


-1 

H 


superhorizon scales 








r H = const.^-- — 










subhorizon scales 




a 



"in a x a \-< fl r-m ^0 

FIG. 3: Evolution of a quantum particle's physical wave 
length A p h = a(t)X as the universe expands. Wave lengths 
which are now of cosmological size were originally minuscule. 
First horizon crossing occurs (at a x ) when A p h becomes com- 
parable to the inflationary Hubble radius ~ c/Hi. At first 
crossing massless, minimally coupled scalars and gravitons of 
that wave length are ripped out of the vacuum by the ex- 
pansion of spacetime. These particles ride the subsequent 
evolution of the universe relatively undisturbed until the sec- 
ond horizon crossing at A p h ~ c/H(t). Then the particles 
manifest themselves as cosmological-scale correlations which 
could not have formed causally after inflation. 

As indiated in Fig. [3] after the first horizon crossing 
particles evolve on superhorizon scales during inflation 
and subsequent eras, to cross an horizon again in the ra- 
diation or matter era. For later use we shall now show 
how to relate the scale factor at the second horizon cross- 
ing to the coordinate wave vector ||fe||. To this end, 
it is convenient to use the redshift, z = ao/a(t) — 1, 
rather than time to label events. This choice implies 
a = an/(l + z). For simplicity, we assume perfect matter 
domination from matter-radiation equality (z cq ~ 3200) 
to the present (z — 0). From Eq. (jTo)l of Sec. [I] we can 
solve for the time in terms of the redshift and express the 
instantaneous Hubble parameter H(t) = 2/Zt in terms of 
z. If we multiply by the scale factor, we find 



aH = ao-ffo Vl + z. (matter domination) (64) 



cq we assume perfect radiation domination. By 



For z > z, 

following the same procedure with Eq. i|15|) . we find 
l + z 



aH = a^Ho 



(radiation domination) (65) 



Suppose galaxies form at z = 10. A physical wave length 
of 10 kpc would have experienced second horizon crossing 
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during the radiation epoch at about Z2 X — 1.5 x 10 7 . 
Therefore the scales relevant to galaxy formation crossed 
during the radiation era, and we conclude that 



\k\ 



a>2 X 



Z2x 



(66) 



We have so far discussed only the fate of virtual par- 
ticles that happen to emerge from the vacuum. The 
amount of particle production also depends upon the rate 
at which they emerge. Most types of massless particles 
possess a symmetry that makes this rate decrease as the 
universe expands. So any of these particles which emerge 
with c||fc p h|| Hi can persist forever without violating 
the uncertainty principle, but not many emerge. 

This symmetry is called conformal invariancc. Particle 
physicists find it convenient to describe symmetries, and 
all the other properties of a dynamical system, in terms 
of the system's Lagrangian. In the classical mechanics of 
point particles the Lagrangian is the difference between 
the kinetic and potential energies, but the concept can 
be generalized to describe any sort of system, including 
the quantum field theories responsible for the vacuum 
polarization effect described in Sec. [H] A field theory is 
said to be conformally invariant if its Lagrangian density 
(the Lagrangian is the spatial integral of the Lagrangian 
density) is unchanged when we multiply each field by a 
certain power of an arbitrary function of space and time 
Q(x). Some interesting fields are the metric g llUl the 
vector potential A^, the Dirac field ipb (with 6 = 1,2,3,4) 
of spin i particles, and the scalar field <j>. Their conformal 
transformations are, 



9^ 



Q 2 9fiv, A^ 



Au, tpb -* H 2 "0b, 



(67) 

A typical conformally invariant Lagrangian density is 
that of electromagnetism, 



Cem = --F a pF pa g a PgP°y 



(68) 



where g = det(g^) denotes the determinant of the metric 
and g^ is its matrix inverse. 

Conformal invariance is so important because there is 
a coordinate system in which a general homogeneous and 
isotropic metric, Eq. (JjjJ, is just a conformal factor times 
the metric of flat space. The change of variables is defined 
by the differential relation, dr\ = dt/a(t), 

ds 2 = -dt 2 + a 2 (t)dx ■ dx = a 2 [-dr] 2 + dx ■ dx] . (69) 

In the (77, x) coordinate system — which we shall hence- 
forth employ, except as noted — the metric and its in- 
verse are, 



g^u 



g^ = a" V 



(70) 



where 



<1 



diag(— 1, 1,1,1). In this coordinate 



same as in flat space when expressed in terms of the con- 
formally rescaled fields Eq. I|67|) with £1 = a -1 . For exam- 
ple, the Lagrangian densities of electromagnetism, mass- 
less Dirac fermions, and a massless, conformally coupled, 
complex scalar are, 



Cem = -~F afi F pa TfP^ , 



(71) 

(72) 
(73) 



where 7^. (with 6, c = 1,2,3,4) are the gamma matrices 
of Dirac theory [14| . 

Although the formalism of Lagrangian quantum field 
theory is a straightforward generalization of quantum 
mechanics, the reader who is unfamiliar with it need not 
fret over how Lagrangians and Lagrangian densities en- 
capsulate the dynamics of any particular system; it is suf- 
ficient for our discussion to note that they do. Therefore, 
any two systems whose Lagrangians agree have the same 
dynamics. We shall exploit this fact at several points 
in the subsequent discussion, starting now. Because the 
massless particle Lagrangian densities in Eqs. l(7T)l - l(75|) 
agree for the conformally rescaled fields with those of flat 
space, it follows that the rate at which virtual particles 
emerge from the vacuum is the same in conformal coor- 
dinates as it is in flat space. Let us call this constant rate 
F. It gives the number of virtual particles emerging per 
unit conformal time r\. This means that the number per 
unit physical time t is, 



dn 
~dt 



drj dn 
dt drj 



(74) 



Hence the emergence rate falls like l/a(i) in physical 
time, and we see that the production of massless, con- 
formally invariant particles is highly suppressed. 

Readers who have studied undergraduate quantum me- 
chanics are familiar with the minimal coupling prescrip- 
tion by which one generalizes the Schrodinger equation 
to account for the interaction of a point particle with 
an electromagnetic field. What one does is to replace 
the derivative operator everywhere with the covariant 
derivative d p + i(e/hc)A fJ ,(x), where e is the particle's 
charge and A p {x) is the 4- vector potential evaluated at 
the particle's position in spacetime. A very similar pro- 
cedure exists for generalizing the dynamics of any sys- 
tem from flat space to curved space. This prescription is 
also called minimal coupling. For electromagnetism and 
massless Dirac fermions it happens to produce the confor- 
mally invariant Lagrangian densities Eqs. (|71() and 17211 . 
However, when a massless scalar is minimally coupled, 
the Lagrangian density which results is not Eq. (|73|l . but 
rather, 



C 



MS 



= -d u 



-cPdupdvW". (75) 



system, a conformally invariant Lagrangian density is the 



This Lagrangian density is not conformally invariant, 
which means that the rate at which massless, minimally 
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coupled scalars emerge from the vacuum per unit physical 
time need not fall off like the rates H74|) of other massless 
particles. We shall now show, by direct calculation, that 
it does not. 

Recall that the Lagrangian of any field theory is ob- 
tained by integrating the Lagrangian density over space. 
Doing this in the original (t, x) coordinate system for 
Eq. I|75|) and applying Parseval's theorem gives 



d 3 xC M s 



(76) 



a 3 (t) 



\\ckl 
« 2 (i) 



m,k)\ 2 }- 



Now note that integration is a form of summation. The 
correspondence can be made explicit by changing to di- 
mensionless variables k — 2irriHi/c, and then exploiting 
the Maclaurin relation between sums and integrals (fa- 
miliar from undergraduate statistical mechanics), 



d 3 k 

(27T) 5 



d 3 n 



E 



(77) 



We can define point particle positions using the real and 
imaginary parts of cf>(t, k), 



<Kt,k) = 



I he 3 

2H] 



(x(t, 



n) + iy(t, n) 



(78) 



Hence the scalar Lagrangian, 
nH ia 3 ^ f - 2 , -2 

Lms^-^-1^{x +y - 



c||fc|| 



\x 2 +y 2 )}, (79) 



reveals that each wave vector of the scalar corresponds to 
two independent harmonic oscillators with the following 
time dependent mass and frequency, 

m(t) =c~ 2 nHia 3 (t) and u(t) = c\\k\\/a(t) . (80) 

Harmonic oscillators with a time dependent mass and 
frequency have been much studied in quantum mechan- 
ics [l9|| . The minimum energy at time t is well known to 
be ^tiu)(t), however, the state with this energy does not 
generally evolve onto itself. For the inflationary case of 
a(t) oc e Hlt , the system's time dependence can be solved 
exactly. The state whose energy is minimum in the dis- 
tant past has instantaneous average (zero-point) energy, 



Eo-^point(t, k) 



hc\\k\\ hH]a(t) 
a(t) 2c\\k\\ 



(81) 



The second term in Eq. (|81l) is attributable to particle 
production. The energy of a single particle with this wave 
vector is 7jc||fc||/a(t), so the average number of particles 
with wave vector k is, 



(82) 



As we expect, N(t, k) is small for very early times and 
becomes comparable to one at horizon crossing. If we 
sum the contributions from all wave vectors that have 
experienced horizon crossing and divide by the spatial 
volume, we find the number density, 



N 
V 



47T 2 C 3 



(83) 



This corresponds to 1/87T 2 particles per Hubble volume 
for each degree of freedom. 

We close by commenting that there can be no ques- 
tion about the reality of inflationary particle production 
because its impact has been detected. There is strong 
evidence that it is what caused the anisotropies imaged 
by WMAP 0. Indeed, all the cosmological structures of 
the current universe are the result of gravitational col- 
lapse into these (originally) quantum fluctuations over 
the course of 13.7 billion years! 



IV. VACUUM POLARIZATION IN INFLATION 

The inflationary Hubble parameter Eq. (1221) corre- 
sponds to an enormous energy, 



hHi ~ 10 3 J ~ 10 13 GeV . 



(84) 



On this scale all the charged particles in the Standard 
Model of particle physics are effectively massless. Even a 
particle we normally consider very massive, such as the 
t quark, has less than 10 -10 times as much rest mass en- 
ergy. However, all but one of the Standard Model charged 
particles are described by the Dirac field ip a , whose La- 
grangian density, Eq. I|72|) . becomes conformally invari- 
ant when we ignore masses. As explained in Sec. IIIII 
conformally invariant particles are not produced much 
during inflation. This means that they do not contribute 
much more to the polarization of the vacuum during in- 
flation than they do in flat space. 

The lone exception within the Standard Model of ele- 
mentary particle theory is the charged sector of the Higgs 
scalar. At low energy it manifests as the longitudinal 
component of the W /± . Its mass of about 80 GeV is also 
insignificant on the scale of inflation. No one really knows 
how it couples to the metric, but the usual assumption, 
based on how the field renormalizes, is minimal coupling. 
We can therefore model it with the Lagrangian density 
of a massless, charged and minimally coupled scalar, 

Csqed = -(d^-ie'A^idv+ie'AJcPg^^ 

= -a 2 (d„ - ie!A tl )4>*{d v + ie'A v )(j)rf v . (85) 



(Here e' = e/hc, and e ~ — 0.30v?ic is the charge of the 
electron. The subscript SQED stands for scalar quan- 
tum electrodynamics.) There may be more so-far undis- 
covered charged scalars of this type lurking between the 
~ 10 2 GeV energies which can be explored at accelerators 
and the enormous energy in Eq. I|84|) of the inflationary 
Hubble parameter. 
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With Ola Tornkvist we have computed the vacuum 
polarization from Csqed Hi! El- ln the x M = (77, x) 
coordinates Eq. (|69[) our result for the polarization bi- 
tensor takes the same form Eq. I|42[) as it does for the 
linear, isotropic medium discussed in Sec. |nj With the 
scale factor normalized to unity at the start of inflation 
the two scalar functions are, 

Xe (x,x') = ^\x,x')-^hx{a)5\x-x') 

D7T 

+ 9 2 |#(A?7)0(At 2 )[1 + ln(if 2 Ar 2 )]} (86) 

tn 2 (x,x') = aa ^J f e(A V )d(Ar 2 )[2 + ln(if f 2 Ar 2 )] 



47T 2 C 4 



(87) 



Here Arj = r\ — rj and At 2 = (77 — i]') 2 — c~ 2 \\x — a?' 1 1 2 . 
xf at (x,x') is the flat space result (|59|) . with i and £' re- 
placed by rj and 77'. This term is renormalized precisely 
as in flat space, and contains no scale factors. The in- 
flationary corrections are completely finite and depend 
on the scale factors a = a(r]) and a' = a(rj'). These 
inflationary corrections come from the long wave length 
virtual particles that are ripped out of the vacuum by the 
inflationary Hubble flow. This should obviously increase 
polarization because it fills spacetime with a plasma of 
charged particles. 

A significant feature of our result is nonzero Sn 2 (x, x') . 
Recall that it must always vanish in flat space quantum 
field theory by virtue of the Lorentz symmetry between 
space and time. The time-dependent metric of inflation, 
Eq. i(7T)|) . does not possess this symmetry, so 5n 2 (x, x 1 ) ^ 
0. In terms of electrodynamics, this means that Csqed 
induces a relative permittivity which is not the inverse of 
the permeability, so the index of refraction is not unity 
even in "empty" space. 

Because the inflationary metric is time dependent, we 
cannot calculate the mass of the photon by checking for 
a pole in the Fourier transform of the susceptibility as 
we did in flat space, Eq. I|48(l . A better way to proceed is 
by comparison with the Proca Lagrangian density which 
governs the dynamics of a fundamental massive photon, 



£ P = —F a0 F pa g a Pg^^—g- 



—^p,A v g v g 



— /V a p ri ap n 

— ^ r ap r pe'l '/ 



27? 



2% 

2 2 

"''<< " - ^-,1,.!., -/''''a-. 



The field equations associated with this Lagrangian den- 
sity are, 

if v ri pa d p F av - c 2 h- 2 m 2 ^" 'A v a 2 = . (89) 

The mass term is distinguished by its factor of a 2 . 

Now recall Maxwell's equations with vacuum polariza- 
tion, Eq. I|41|l . which we rewrite without the current, 



^ v T] pa d p F av {x) + / dVpiHs, x')A v (x') = 0. (90) 



We also recall the polarization bi-tensor Eq. 142|) . 

TTl v tx,x') = -[^ v r,P a -^Pr," v ]d' p d a xe{x,x') 

+ r 1 » i rf j [5 ij d' k d k -d' i d j ]5n 2 {x,x'). (91) 

We see from Eq. (I87|l that 5n 2 (x, x') contributes a factor 
of a 2 . The Xe{x,x') term, Eq. (jHEJ), has at most a single 
factor of a, but note from d a a = 5®Hia 2 that this term 
can also give an a 2 in Eq. I|9U|) . A comparison with the 



Proca equations (|8T)|) suggests 



y/aHih/c 2 . 



We can obtain a quantitative result by solving Eq. l|9"U|) 
perturbatively in a. First expand the vector potential in 
a series of terms A^\x) which go like a", 



A,(x)=A^\x)+AW(x) + 



(92) 



Now recall that the polarization bi-tensor is first order 
in a, and it groups the terms in Eq. (|90|l in powers of a. 
We see that Ap°\x) obeys the classical equation, 



rrd p F<$(x) - 0, 



(93) 



the general solution of which consists, in the Coulomb 
gauge, of a superposition of transverse plane waves, 

A { ° ] (x) = e^e- 1 ^'^ 3 , where e - = he,. 

(94) 

The order a correction obeys, 
t v r } P"d p F^{x)+ [ d 4 x'ru-](x,x')Ai°\x') = 0. (95) 



Now substitute Eq. (|94|l and evaluate the integral as- 
suming the photon experienced first horizon crossing (see 
Fig. |2| long before, and after a long period of inflation, 



c||fc|| , 



(96) 



After some tedious expansions, the result is, 

r d i x'[^Il v }{x,x l )A^\x l ) (97) 



c\\k\ 



-lnl 

7T V Hr 



+ 0(1) a 2 + 0(a). 



The analogous first order Proca equation, 

V ^ V pa d p F^ - c^m^A^a 2 = 0, 
implies that the photon mass must be, 



\fac 2 hHj 



H«)+o (1 , 



(99) 



We have thus discovered that the photon of scalar elec- 
trodynamics acquires mass during inflation. 
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V. HARTREE APPROXIMATION 

A simple way of obtaining almost the same result was 
previously suggested by one of us (Prokopec) in col- 
laboration with Anne-Christine Davis, Konstantinos Di- 

mopoulos, and Ola Tornkvist HHHI^- The technic l ue 
is to pretend that photons move in the quantum me- 
chanical average of the scalar field. This is known as the 
Hartree or mean field approximation. First year graduate 
students should be familiar with its use in quantum me- 
chanics to treat interacting, multi-electron atoms, and in 
statistical mechanics to treat interacting, multi-particle 
systems such as gases subject to the van der Waals force. 

To implement the Hartree approximation we first take 
the average of Lsqed in Eq. I|85|l over quantum mechani- 
cal fluctuations of the scalar field. Of course this average 
eliminates the scalar fields, but it leaves behind some 
function of the vector potential, 



-SQED 

- ^(ff^fTa 3 . (100) 

Now add this function to the Maxwell Lagrangian density 
Eq. PI). 

By using the sophisticated techniques of quantum field 
theory, we can show that the quantum average of the 
scalar's norm-squared consists of a divergent constant 
plus a finite term that grows like the logarithm of the 
scale factor Hal: 



4>* (x)(j>(x) 
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(101) 



The other averages in Eq. I|100|l are either zero or else 
they do not multiply functions of the vector potential. 
The Hartree approximation Lagrangian density is there- 
fore, 
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U.V. + — T — In(o) A^A^a 2 . (102) 



Expression (|102l) contains ultraviolet divergences be- 
cause virtual particles of arbitrarily large wave vector 
contribute to the average. This is the same origin as the 
divergences we found in Sec.[H] and it is the ultimate ori- 
gin of all ultraviolet divergences in quantum field theory. 
The only thing that need concern us here is the depen- 
dence of the divergent terms on the electromagnetic vec- 
tor potential A^x). The divergence without any vector 
potentials is harmless, but the other one could only be 
renormalized using a fundamental photon mass, which 
we do not have. This is one reason why the vacuum 
polarization — which can be consistently renormalized 



- is the correct way to study the kinematical proper- 
ties of photons. But let us simply ignore the divergences 
in Eq. p()2|l . A comparison of the finite parts with the 
Proca Lagrangian density Eq. suggests the corre- 
spondence, 



\J~OLC 2 TiHj 
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(103) 



Complete agreement with Eq. 1|99|) requires only the addi- 
tional assumption that the growth of Eq. (|103fl ceases for 
the mode of wave vector k (which we assume obeys || k || ^> 
Hi/c) when it experiences horizon crossing, a — c\\k\\/Hi 
(see Fig. |3J|. This point of view is consistent with the 
causal picture according to which a photon's mass only 
receives contributions from virtual scalars whose wave 
lengths are greater than the photon's wave length. 

We conclude this section by commenting on the size 
of the photon mass induced by our mechanism. During 
inflation we have m 7 ~ 10 13 GeV/c 2 , which is enormous 
compared to the center-of-mass energies of ~ 100 GeV 
attainable in the largest accelerators. A photon mass is 
not detected today because our result derives from the 
huge density of free charged particles ripped out of the 
vacuum by inflation. This plasma has been thoroughly 
dissipated at any wave length we can access in today's 
laboratories. 

According to the supernovae results 0, the cur- 
rent universe may be entering another phase of inflation. 
This late inflationary phase also leads to a nonzero pho- 
ton mass, but with the replacement of Hj by the vastly 
smaller Hubble parameter of today, Ho. This substitu- 
tion in Eq. (|99|l gives a minuscule photon mass, m 7 ~ 
10~ 42 GeV/c 2 , which is far below the best current labora- 



tory bounds of m 7 < 10" 49 kg w 10" 23 GeV/c 2 [ljJll 



VI. COSMOLOGICAL MAGNETIC FIELDS 

The phenomenon of nonzero photon mass during in- 
flation offers a fascinating 4-dimensional analogue to the 
Schwinger model in two dimensions |l7j| . However, it was 
proposed [22,|2i|,|3] not for aesthetic appeal, but rather 
to explain the curious fact that galaxies seem to possess 
magnetic fields, which are correlated on scales of a few 
kilo-parsecs, and whose strength is typically a few micro- 
Gauss 26]. (The conversion factors to MKSA units are 
1 J/m 3 = 10 Gauss 2 and lOkpc ~ 3.1 x 10 20 m.) There is 
also evidence that galactic clusters possess micro-Gauss 
magnetic fields correlated on scales of 10-100 kpc |27j . 
Although some of the material in this section is more 
technical than before, wc present it to illustrate how the 
phenomenon of vacuum polarization during inflation may 
have left an observable consequence. 

The difficulty with cosmic magnetic fields is not their 
field strengths but rather their enormous coherence 
lengths. A galaxy's differential rotation can combine 
with the turbulent motion of ionized gas to power a phe- 
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nomenon known as the a-ui dynamo 28]. In this mech- 
anism the lines of a coherent seed field are stretched by 
rotation, twisted by turbulence, and then recombined 
to result in an exponential amplification, B(t) oc e*/ r . 
Kinetic energy from the turbulent motion is converted 
into magnetic field energy in this way until equiparti- 
tion is reached. Although many astrophysicists question 
the a-uj dynamo, it is significant that the measured field 
strengths are at roughly the equipartition limit [22} ■ 

There is no general agreement on a reasonable value 
for the dynamo time constant t; estimates vary from 
0.2 to 0.8 billion years |2^. By observing a surpris- 
ingly large polarization in the cosmic microwave back- 
ground photons, the WMAP satellite has seen reioniza- 
tion from the first star formation at about 0.2 billion 
years into the 13.7 billion years of the universe's exis- 
tence Q . One might expect large spiral galaxies to form 
at about 0.4 billion years [2^|, which implies dynamo op- 
eration for 13.7 — 0.4 = 13.3 billion years, or between 17 
to 66 time constants. Exponentiation results in amplifi- 
cation factors ranging from e 17 ~ 2.4 x 10 7 to as large 
as e 66 ~ 4.6 x 10 28 . Therefore the cosmological magnetic 
fields of today might derive from correlated seeds as weak 
as 10~ 34 Gauss at the time of galaxy formation. The real 
question is what produced the correlated seed fields in 
the hot, dense and very smooth early universe? 

In the following, we argue that the nonzero photon 
mass of inflation might help to answer this question. As 
explained in Sec. II I II a nonzero mass suppresses the cre- 
ation of particles. On the other hand, it vastly enhances 
the zero-point energy that quantum mechanics predicts 
must reside in each photon wave vector k, even if there 
are no particles with that wave vector anywhere in the 
universe. With no mass this zero-point energy falls as 
the universe expands, 



q and p stand for the position and momentum operators 
of this oscillator. The Hamiltonians before and after are, 

Hb = ^ + JmfiV , H A = %- + \mo?q 2 . (106) 
Am 2 Ira A 

Before transition the system is in its ground state, 



H B \0) = -hQ\0) 



(107) 



This result can be obtained from the standard relations, 

q = \^n(«-^)> p = ^(-- + - t ),d08) 

where a — aoexp(— i£lt) and = ajexp(?f2t) are the 
lowering and raising operators, respectively, with ao|0) = 
a|0) = 0, not to be confused with the scale factor a(t). 
The kinetic and potential terms each contribute half, 
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After the transition the system is no longer an eigenstate, 
but we can find its average energy from the fact that the 
expectation values of p 2 and q 2 are continuous, 
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(110) 



If we compare this energy with the post-transition eigen- 
states (N + ^)nu>, we see that the average occupation 
number after transition is, 
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A nonzero photon mass causes the zero-point energy of 
wave vectors that have experienced first horizon crossing 
to approach a constant instead, 



E^(t, k) = - 
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(105) 



After the end of inflation this wave vector eventually ex- 
periences second horizon crossing, c||fc|| ~ a(t)H(t). If 
the mass goes to zero quickly thereafter, about half of the 
enormous zero-point energy must be shed in the form of 
long wave length photons at numbers vastly higher than 
thermal. The idea is that the mysterious seed fields de- 
rive from these long wave length photons becoming frozen 
in the plasma of the early universe. 

Consider a wave vector k that is about to experi- 
ence second horizon crossing. Each polarization of this 
system behaves as an independent harmonic oscillator 
whose frequency is suddenly changed from a large value 
il = m-fC 2 /Ti to a much smaller one u> = c\\k\\/a(t). Let 



The substitution in Eq. (|lll|l of the scale factor at the 
second horizon crossing, a = ci2x as given in Eq. 16611 . 
results in the average occupation number for each polar- 
ization of wave vector k, 



N(k) 



Jz^h\\k\\'< 



(112) 



(Recall that k p h — ||fc||/ao is the physical wave num- 
ber measured today.) To within factors of order unity, 
the temperature at time t is T — Toao/a(t), where 
To ~ 2.73 K is the current temperature of the cosmic mi- 
crowave background. The number of thermal photons of 
wave number k and a fixed polarization obeys the Planck 
distribution, 



N th (k) 



exp 



\kRT a J 



ir fceTpao 



(113) 



where ks — 1-38 x 10 - 23 J/K is the Boltzmann constant 
and the final relation on the right applies in the long wave 
length (IR) limit. The ratio of m 7 photons to thermal 
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ones can be expressed in terms of the present-day physical 
wave length A , 



N(k) 
N th (k) 
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(114) 



Recall that we normalize the scale factor to one at the 
start of inflation. For models with a long period of in- 
flation the final factor in square brackets is dominated 
by ln(ao), which might be quite large. We parameterize 
our ignorance of the number of inflationary e-foldings by 
defining, 
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r(10 kpc), 

If we work out the other numbers, we obtain, 



N(k) 
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(115) 



(116) 



We see that m 7 photons are negligible compared to ther- 
mal ones on the Ao ~ 0.002 m scale of the cosmic mi- 
crowave background, but they are enormously dominant 
on the Ao ~ lOkpc scale relevant to galaxies. 

At first the energy in these photons is almost com- 
pletely electric, but Maxwell's equations carry it to the 
magnetic sector. The physical magnetic field in a homo- 
geneous and isotropic geometry is, 



B\t,x) 



e^ k F jk (t,x)/a 2 (t) 



(117) 



If we assume that half of the energy of the massive pho- 
tons winds up in these magnetic fields, we conclude that 
their spatial Fourier transforms obey, 



W ( / . k)W ( / . q) ) = (2,r) 3 * 3 (£+gW£)M*l[, (118) 
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The quantity of interest is the magnetic field averaged 
over a region of coordinate size £ = £o/a(t), 

B%x;l ) = {2ir£ 2 yi (<$ye-^B\t,y), (120) 
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B l (t,x;£o) is an operator, but its average is a number, 



B 2 (t,£ ) = (B l (t,x;£ Q )B\t,x-£ )) 
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If we plug in the known numbers, we find 



B(t, £ Q )~^fex 10~ 33 Gauss 



l + z 
4/10 kpc 



(122) 



(123) 



This value is already within the lower range of conceiv- 
able seed fields. Turbulent evolution might contribute a 
factor of ten by transferring power from small scales |22| . 
An additional factor of (p ga i/po) 3 — 3 x 10 3 accrues 
from field compression when the proto-galaxy collapses. 
If we assume e ~ 1 and galaxy formation at z ~ 10, 
we might expect field strengths of about 10~ 28 Gauss at 
£ - lOkpc. 

It should be emphasized that this is just one of 
many potential explanations for cosmological magnetic 
fields |27|. This section's analysis is also highly simpli- 
fied. We need to better understand the process through 
which a given wave vector's mass dissipates at second 
horizon crossing. A proper calculation would also re- 
quire careful study of the dynamics of the electric and 
magnetic fields during the epochs of reheating, radiation 
domination, and matter domination. 
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